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The CLEO collaboration has studied two-body charmless hadronic decays of B mesons into final 
states containing two pseudo-scalar mesons, or a pseudo-scalar and a vector meson. We summarize 
and discuss results presented during the winter/spring 1999 conference season, and provide a brief 
outlook towards future attractions to come. 

In particular, CLEO presented preliminary results on the decays B — > p° (Br(B ± — > 7r ± p°) = 
(1.5±0.5±0.4)xl(T 5 ), B -> 71-=^ (Br(B -> tt^) = (3.5±i ±0.5) x 1(T 5 ), B -> n ± K* :f (Br(B -> 
n ± K* T ) = (2.2±° %t° A 5 ) x HT 5 ), and B ± -» K^ (Br(B ± -> iTV) = (1.5 ± 0.4 ± 0.3) x 10~ 5 ) 
at DPF99, APS99, APS99, and ICHEP98 respectively. None of these decays had been observed 
f*"^ ' previously. The first two of these constitute the first observation of hadronic b — > u transitions. In 

addition, CLEO presented preliminary updates on a large number of previously published branching 
■ fractions and upper limits. 

0\ 
On 

X • I. INTRODUCTION AND OVERVIEW 

Oh. The phenomenon of CP violation, so far observed only in the neutral kaon system, can be accommodated by 
a complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark- mixing matrix jjj. Whether this phase is the 
correct, or only, source of CP violation awaits experimental confirmation. B meson decays, in particular charmless 
B meson decays, will play an important role in verifying this picture. 

The decays B — > iy + tt~ and B — > p + ir~ , dominated by the b — > u tree diagram (Fig. 0(a)), can be used to measure 
CP violation due the interference between B° — B° mixing and decay. However, theoretical uncertainties due to the 
presence of the b — > dg penguin diagram (Fig. |l](b)) ("Penguin Pollution") make it difficult to extract the angle a of 
the unitarity triangle from B — > tt + tt~ alone. Additional measurements of B — > tt^tt , B,B — > 7r°7r°, or a flavor 
tagged proper-time dependent full Dalitz plot fit for B — > 7r + 7r~7r° and the use of isospin symmetry may resolve 
these uncertainties || || |Q. Alternatively, measurement of CP violation due to the interference of B —> 7r ± XcO 
and B — > ir^p in B ± — * ir+n^^ || may provide information about the angle 7 of the unitarity triangle. Neither 
flavor tagging nor a measurement of the proper-time before decay of the B meson is required in this case. Extraction 
of the angle 7 from this measurement is subject to theoretical uncertainties due to Penguin Pollution. However, 
factorization predicts this to be less severe here than in B — > 7r+7r~ due to at least partial cancelation of the gluonic 
penguin contribution among short distance operators with different chirality. 

B — * Kit decays are dominated by the b — ► sg gluonic penguin diagram, with additional contributions from b — > u 
tree and color-allowed electroweak penguin (Fig. |l](d)) processes. Interference between the penguin (Fig. |l|(b),(d)) and 
spectator (Fig. |l](a) ,(c)) amplitudes can lead to direct CP violation, which would manifest itself as a rate asymmetry for 
decays of B and B mesons. Several methods of measuring the angle 7 using only decay rates of B — > Kir, tttt processes 
were also proposed ||. This is particularly important, as 7 is the least known parameter of the unitarity triangle and 
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is likely to remain the most difficult to determine experimentally. The ratios R = B(B — > K ± n T )/B(B ± — > if ?? 1 ' 1 ) 
0, and i?* = B(B ± — > K a ir ± )/2B(B ± — > K^ir ) ||], were recently suggested as a way to constrain 7. Electroweak 
penguins and final state interactions (FSI) in _B — > i^7r decays can significantly affect the former method || , whereas 
the latter method requires knowledge of the ratio \ (T+C)/ P\ s of spectator to penguin amplitudes in b — ► s transitions. 
Uncertainties due to FSI and electroweak penguins are eliminated using isospin and fierz-equivalence of certain short 
distance operators. Studies of B decays to KK final states can provide useful limits on FSI effects fTofl . 
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FIG. 1. The dominant decay processes are expected to be (a) external W-emission, (b) gluonic penguin, (c) internal 
W-emission, (d) external electroweak penguin. 

B decays to r]'K®, p°K®, and 4>K° may allow for future measurements of sin 2/3, (3 being the third angle of the 
unitarity triangle. This is of interest because one probes the interference between the amplitudes for b — » s penguin 
and B — B° mixing, rather than b — > c tree and B° — B° mixing, as done in the more ubiquitous B — » ipK® decay. It 
has been argued jllj that a variety of new physics scenarios would affect the CP violating phase of the b — > s penguin 
only, leaving the phases of B° — B° mixing and b — > u tree amplitudes unchanged. Such new physics scenarios would 
thus lead to a difference between proper-time dependent CP violation as measured for example in B decays to rj'K® 
as compared to B — > ipK® . 

The present paper presents preliminary CLEO results on two-body charmless hadronic decays of B mesons into 
final states containing two pseudo-scalar mesons (B — * PP), or a pseudo-scalar and a vector meson (B — > PV). 
Section II discusses the analysis technique that is common to all of these analyses. Results on B — > PP and B — > PV 
are presented in Sections III. Section IV discusses possible implications of some of the measurements presented. 

II. DATA ANALYSIS TECHNIQUE 

The data set used in this analysis is collected with the CLEO II and CLEO II. 5 detectors at the Cornell Electron 
Storage Ring (CESR). Roughly 2/3 of the data is taken at the T(4S) (on-resonance) while the remaining 1/3 is taken 
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just below BB threshold. The below-threshold sample is used for continuum background studies. The on-resonance 
sample contains 5.8 million BB pairs for all final states except p + h~ , K* + h~ (h + being a charged kaon or pion), 
and p°K®. For those final states a total of 7.0 million BB pairs was used. This is an 80% increase in the number 
of BB pairs over the published analyses In addition, we have re-analyzed the CLEO II data set with improved 
calibration constants and track-fitting algorithm allowing us to extend our geometric acceptance and track quality 
requirements. This has lead to an overall increase in reconstruction efficiency of 10 — 20 % as compared to the 
previously published analyses. The CLEO detector has been decommissioned for a major detector and accelerator 
upgrade. Preliminary results based on the full data set of roughly 10 million BB pairs are expected to be ready for 
the summer conferences in 1999. 

CLEO II and CLEO II. 5 are general purpose solenoidal magnet detectors, described in detail elsewhere Q. In 
CLEO II, the momenta of charged particles are measured in a tracking system consisting of a 6-layer straw tube 
chamber, a 10-layer precision drift chamber, and a 51-layer main drift chamber, all operating inside a 1.5 T super- 
conducting solenoid. The main drift chamber also provides a measurement of the specific ionization loss, dE/dx, 
used for particle identification. For CLEO II. 5 the 6-layer straw tube chamber was replaced by a 3-layer double sided 
silicon vertex detector, and the gas in the main drift chamber was changed from an argon-ethane to a helium-propane 
mixture. Photons are detected using a 7800-crystal CsI(Tl) electromagnetic calorimeter. Muons are identified using 
proportional counters placed at various depths in the steel return yoke of the magnet. 

Charged tracks are required to pass track quality cuts based on the average hit residual and the impact parameters 
in both the r — (f> an d r — z planes. Candidate K s are selected from pairs of tracks forming well measured displaced 
vertices. Furthermore, we require the K s momentum vector to point back to the beam spot and the ir + ir~ invariant 
mass to be within 10 MeV, two standard deviations (a), of the K s mass. Isolated showers with energies greater than 
40 MeV in the central region of the Csl calorimeter and greater than 50 MeV elsewhere, are defined to be photons. 
Pairs of photons with an invariant mass within 25 MeV (~ 2.5cr) of the nominal ir° mass are kinematically fitted with 
the mass constrained to the ir° mass. To reduce combinatoric backgrounds we require the lateral shapes of the showers 
to be consistent with those from photons. To suppress further low energy showers from charged particle interactions 
in the calorimeter we apply a shower energy dependent isolation cut. 

Charged particles are identified as kaons or pions using dE/dx. Electrons are rejected based on dE/dx and the 
ratio of the track momentum to the associated shower energy in the Csl calorimeter. We reject muons by requiring 
that the tracks do not penetrate the steel absorber to a depth greater than seven nuclear interaction lengths. We 
have studied the dE/dx separation between kaons and pions for momenta p ~ 2.6 GeV/c in data using £>* + -tagged 
D° -> R-TT+ decays; we find a separation of (1.7 ± 0.1) a for CLEO II and (2.0 ± 0.1) a for CLEO II.5. 

We calculate a beam-constrained B mass M = \J E^ — p 2 B , where ps is the B candidate momentum and Eh is the 
beam energy. The resolution in M ranges from 2.5 to 3.0 MeV/c 2 , where the larger resolution corresponds to decay 
modes with a high momentum ir°. We define AE — J2i Ei — Eh, where Ei are the energies of the daughters of the B 
meson candidate. The resolution on AE is mode-dependent. For final states without 7r°'s the AE resolution for CLEO 
11(11.5) is ~ 20— 26(17— 22)MeV. For final states with a high momentum 7r° the AE resolution is worse approximately 
by a factor of two and becomes asymmetric because of energy loss out of the back of the Csl crystals. The energy 
constraint also helps to distinguish between modes of the same topology. For example, AE for B — > K + tt~ , calculated 
assuming B — > tt + it~ , has a distribution that is centered at —42 MeV, giving a separation of 1.6(1.9)<7 between 
B — ► K + ir~ and B — > tt + it~ for CLEO 11(11.5). In addition, AE is very powerful in distinguishing B — > K*°ir + from 
B — > p°ir + , especially if the positive track from the vector meson is of low momentum. 

We accept events with M within 5.2 — 5.3 GeV/c 2 . The fiducial region in AE depends on the final state. For 
B — ► PP we use \AE\ < 200(300) MeV for decay modes without (with) a ir° in the final state. The selection criteria 
for B — ► PV are listed in Table II. This fiducial region includes the signal region, and a sideband for background 
determination. 

We have studied backgrounds from b — > c decays and other t->« and b — > s decays and find that all are negligible 
for B decays to two pseudo-scalar mesons. In contrast, some of the B decays to a pseudo-scalar and a vector meson 
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have significant backgrounds from b — > c as well as other charmless B decays. We discuss these in more detail below 
in Section III. However, the main background in all analyses arises from e + e~ — > qq (where q — u,d, s, c). Such events 
typically exhibit a two-jet structure and can produce high momentum back-to-back tracks in the fiducial region. To 
reduce contamination from these events, we calculate the angle 9s between the sphericity axis juj of the candidate 
tracks and showers and the sphericity axis of the rest of the event. The distribution of cos 9s is strongly peaked at 
±1 for qq events and is nearly flat for BB events. We require | cos 9s \ < 0.8 which eliminates 83% of the background 
for all final states except those including rf or <j). For the latter final states a looser cut of | cos 9s\ < 0.9 is used. 

Using a detailed GEANT-based Monte-Carlo simulation JlEj] we determine overall detection efficiencies {£ ) ranging 
from a few % to 53% in B — > K + n~ . Efficiencies are listed for all decay modes in the tables in Section III. We 
estimate systematic errors on the efficiencies using independent data samples. 

Additional discrimination between signal and qq background is provided by a Fisher discriminant technique as 
described in detail in Ref. p6| . The Fisher discriminant is a linear combination T = 53i=i a iVi where the coefficients 
cdj are chosen to maximize the separation between the signal and background Monte-Carlo samples. The 11 inputs, j/j, 
are | cos0 con d| (the cosine of the angle between the candidate sphericity axis and beam axis), the ratio of Fox- Wolfram 



moments H2/HQ |17|, and nine variables that measure the scalar sum of the momenta of tracks and showers from the 
rest of the event in nine angular bins, each of 10°, centered about the candidate's sphericity axis. Some of the analyses 
(final states including rf or (f>) use |cos#s| (the angle between the B meson momentum and beam axis) instead of 
H2I ' Hq as one of the inputs to the Fisher discriminant. 

We perform unbinned maximum- likelihood (ML) fits using AE, M, J 7 , \cos9b\ (if not used as input to J 7 ) and 
dE/dx (where applicable) as input information for each candidate event to determine the signal yields. Resonance 
masses (7/ and vector resonances) and helicity angle of the vector meson are also used as input information in the fit 
where applicable. In each of these fits the likelihood of the event is parameterized by the sum of probabilities for all 
relevant signal and background hypotheses, with relative weights determined by maximizing the likelihood function 
(£). The probability of a particular hypothesis is calculated as a product of the probability density functions (PDFs) 
for each of the input variables. Further details about the likelihood fit can be found in Ref. |l(|. The parameters for the 
PDFs are determined from independent data and high-statistics Monte-Carlo samples. We estimate a systematic error 
on the fitted yield by varying the PDFs used in the fit within their uncertainties. These uncertainties are dominated 
by the limited statistics in the independent data samples we used to determine the PDFs. The systematic errors on 
the measured branching fractions are obtained by adding this fit systematic in quadrature with the systematic error 
on the efficiency. 

In decay modes for which we do not see statistically significant yields, we calculate 90% confidence level (C.L.) 
upper limit yields by integrating the likelihood function 
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J? C m ^(N)dN 
f °° C max (N)dN W 

where C max (N) is the maximum C at fixed N to conservatively account for possible correlations among the free 
parameters in the fit. We then increase upper limit yields by their systematic errors and reduce detection efficiencies 
by their systematic errors to calculate branching fraction upper limits given in Table I and IV. 



III. RESULTS 

Given the enormous number of results to summarize in this Section, we choose to show figures only for those decay 
modes for which we observe statistically significant yields, and no branching fraction measurements have previously 
been published. Additional figures for preliminary updates on previously published branching fraction measurements 
can be found elsewhere. [fl8[ 

The figures we show are contour plots of — 21n£ for the ML fit as well as projection plots for some of the fit 
inputs. The curves in the contour plots represent the ncr contours , which correspond to the increase in — 21n£ 
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by n 2 . Contour plots do not have systematic errors included. The statistical significance of a given signal yield is 
determined by repeating the fit with the signal yield fixed to be zero and recording the change in — 21n£. For the 
projection plots we apply additional cuts on all variables used in the fit except the one displayed. These additional 
cuts suppress backgrounds by an order of magnitude at signal efficiencies of roughly 50%. Overlaid on these plots are 
the projections of the PDFs used in the fit, normalized according to the fit results multiplied by the efficiency of the 
additional cuts. All results shown are preliminary. Not all published analyses [M have been updated yet. 

A. B Decays to Two Pseudo-scalar Mesons 

Table 1 lists the preliminary CLEO results for B decays to two pseudo-scalar mesons. Not all possible final states 
with two pseudo-scalar mesons have been updated yet. For published results please refer to Ref. fl^| . 

Figure || illustrates a contour plot for the ML fit to the signal yield (N) in the track n° final state. The dashed 
curve marks the 3cr contour. To further illustrate the fit, Figure || shows M (AE) projections as defined above. 
Events in Figure || are required to be more likely to be kaons than pions according to dE/dx. We find statistically 
significant signals for the decays B — > A" ± 7r =F , B^ —> K^ir , B ± — > Kgir^, as well as the two B — ► rf ' K decays. The 
corresponding branching fractions are listed in Table |. Table | also shows 90% confidence level upper limits for all 
the decay modes where we do not measure statistically significant yields. 



TABLE I. 


Summary of preliminary CLEO results 


for B decays to two pseudo- 


-scalar mesons. 


Yield and efficiencies in decay 


modes including rj' refer to rj' — * rjir + ir , rj — ► *yy (rj' 


— > pj) decays. 






Mode 


Eff (%) 


Yield 


Signif 


BR/UL (10- b ) 




53 ±5 


43.1 +y 8 ^ 


> 6(7 


1.4 ±0.3 ±0.2 


K + TT 


42 ±4 


38.1 +9 8 '£ 


> 6cr 


1.5 ±0.4 ±0.3 


K ^ 


15 ±2 


12 3 +47 


> 5a 


1.4 ±0.5 ±0.2 


ri'K* 


5(11) 


18.4(50.2) 


> 6a 


7.4±°J ± 1.0 


V 'K° 


1.6(3.4) 


5.4(12.7) 


> 5a 


5.9±i;| ± 0.9 




53 ±5 


11.5j$ 


< 3a 


< 0.84 




42 ±4 


14 q+ 81 


< 3a 


< 1.6 




5(11) 


1.0(0.0) 




< 1.2 


K*K* 


53 ±5 






< 0.23 




15 ±2 






< 0.93 
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FIG. 2. Contour of the -21n£ for the ML fit to N K±7! o and N v ±„o for B ± -> ^tt and -> tt+tt . 




-0.20 -0.10 0.00 0.10 0.20 

AE (GeV) 

FIG. 3. Projection plots in B ± -> /S^tt . AS for B± -> if* 7T° is centered around — 42MeV because we use pion mass when 
calculating the energy of a track. Only events for which the candidate track is more likely to be a kaon than a pion according 
to dE/dx enter these figures. 
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B. B Decays to a Pseudo-scalar and a Vector Meson 



Hclicity conservation dictates that the polarization of the vector in B — > PV is purely longitudinal (helicity = 
state). The kinematics of these decays (assuming two-body decay of the vector) therefore results in a final state with 
two energetic particles and one soft particle. The pseudo-scalar P is always very energetic, with a momentum range 
from 2.3 to 2.8 GeV. On the other hand the decay daughters from the vector meson have a very wide momentum range. 
While the more energetic particle has momentum between 1.0 and 2.8 GeV, the soft particle can have momentum as 
low as 200 MeV. 

The backgrounds from BB events are potentially dangerous as they may peak in either or both of the M and AE 
distributions. There are two types of BB backgrounds that can contribute to PV: b — > c processes and other rare b 
processes. 

Among the B — > PV modes we are searching for, B — > n^V and B — > K ± V can be well separated, using the dE/dx 
information of the very energetic 7r + or K + and the separation in AE, just like the B — ► PP modes. 

Crosstalk of two kinds exist among PV modes. First, p <-+ K* misidcntification is possible for track ir° as well 
as two track decays of the p or K* if the fast particle is misidentified due to the limited particle ID for fast tracks. 
Crosstalk among B + — > h + p° and B + — > h + K*° can be controlled to a level of 20% or less just by requirements on 
dE/dx (2 ct) of the decay daughters of the vector meson. Further separation is achieved by using AE and resonance 
mass of the vector as inputs to the likelihood fit. Second, it is possible to swap a slow momentum pion from the vector 
with a slow momentum pion from the other B. This is particularly severe for slow momentum ir°, as the fake/real 
7T° ratio is about a factor 20 worse for the slow pions than the fast pions from the vector. In such cases we impose 
helicity requirement to remove the region with soft 7r°. Using data (doubly charged vector candidates) and Monte 
Carlo we determine the remaining backgrounds from other rare b processes to be small effects that we correct for. 

The dominant b -> c background for PV{h+h+h-) is B+ -> D°ir+, where D° -> tt+tt" or D° -> K + n~. This 
particular background has exactly the same final state particles as the PV(h + h + h~) signal and therefore peaks both 
in M at 5.28 GeV and in AE at 0.0 GeV. Other b -» c processes B+ -» L»° i o+ and B° -> L»*"7r+ will have peak 
structure in M, but not in AE due to the missing soft particle. Because of the large b — > c branching ratio, the 
contribution from these processes needs to be highly suppressed. We apply a D° (30MeV) veto to all possible h + h~ 
combinations in PV(h + h + h~) modes. 

Similarly the b — > c background for PV(h + h~Tr°) are _B° — > D~7r + where D~ — > 7r~7r or _B° — ► _D°7r° where 
D° — > K~tt + . However, their contribution is negligible due to the branching ratios involved. 

Finally, there are potentially backgrounds from non-resonant B decays to three-body final states. We test for such 
backgrounds in data by allowing a non-resonant signal contribution in the fit, as well as by determining the fit yield 
in bins of helicity angle. Neither of these tests shows any evidence of non-resonant contributions to any of our final 
states. The increase of the error on the fitted yield due to possible non-resonance contributions is accounted for as 
part of our systematic errors. 

1. First Observation of B ± — ► n ± p" 

We select separate h + p° and h + K*° samples as discussed above and in Table II. We then fit for the B + — > ir + p° and 
B + — ► K + p° components in this h + p a sample, as well as a B — ► ir + K* a reflection, averaging over charge conjugate 
modes. Similarly we select a h + K*° sample and fit for the B + — > ir + K* a and B + — > K + K*°, as well as a B + — ► ir + p° 
reflection. We do not attempt a simultaneous fit to the h + p° and h + K* a samples at this point as this would require 
us to model the full momentum dependence of AE, dE/dx, and resonance mass in order to separate p and K* 
contributions. 

The variables M, T, E(7T7r7r) - E b (E^Tr) - E 6 ), dE/dx of h in B -» hp (B -» hK*°), Mass of p° (K*°) 
candidate and cos(p° (K*°) helicity angle) are used to form probability density function (PDF) to perform the ML 
fit for B ± — ► h ± p° (B^ 1 — » li^K*®) sample. We do not use dE/dx for the daughters of the vector meson in the fit. 
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Efficiencies and results are summarized in Table IV. A significant signal in _B ± — > 7r ± j o° is observed. The contour 
and projection plots are shown in Fig. 0. 



Contour Plot B + -» h>° where p" -» -n + -n~ B + -> h>° where p° -» tt+tt 




10 20 30 40 50 60 5.200 5.225 5.250 5.275 5.300 0.57 n - 77 97 

B + ^>7T + p Yield Mbc in GeV p Mass 

FIG. 4. Contour and projection plot for — > ir ± p°. 



TABLE II. Event Selection for B — > PF decays 



Sample 


E — Ebeam 


Resonance Mass Window 


Cos(Resonance Helicity Angle) 


fc=V 


|E(7r7rvT)-Ebcam| <100MeV 


200MeV 


-0.9 - 0.9 


h±K* u 


|E(7rJf7r)-Ebeam| <100MeV 


75MeV 


-0.9 - 0.9 


h ± p+ 


|E(7r7r7r u )-Ebeam| <150MeV 


200MeV 


0.0 - 0.9 




|E( 7 rA'7r u )-Ebeam| <300MeV 


200MeV 


0.1 - 1.0 




E(7rif7r u )-Ebeam| <200MeV 


200MeV 


-0.86 - 1.0 



The contribution of b — > c and other related rare B processes are small but not negligible. They are evaluated 
using about 25 million generic bb Monte Carlo events, and specific Monte Carlo samples for all the rare B processes 
mentioned in this paper. The dominant contributions are listed in Table III. All other contributions are negligible. 



TABLE III. Contributions to the B — ► ir p° Yield from Non-continuum physics backgrounds 



Decay Process 


Contribution to B ± -> i^p" Yield 


b -» c 


0.9±0.7 


B [) - Tv+p- 


0.7±0.3 


B u - K+p- 


O.liO.l 




0.3±0.2 


TOTAL 


2.0±0.8 
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The final B ± — > ^p yield after background subtraction is: 26.1_g events, leading to a branching fraction 
measurement of B(B + — > tt + p°) — (1.5 ± 0.5 ± 0.4) x 10~ 5 . This is the first observed hadronic b — > it transition. 

J2. First Observation of B° — > 7r ± p :F 

As discussed above, the 7r° daughter of the p has a bi-modal momentum distribution due to the longitudinal 
polarization (helicity = 0) of the p. The ratio of real to fake ir° is roughly 1/2 for the low and 10/1 for the high 
momentum 7r° region. This leads to largely increased backgrounds from all sources as well as multiple entries per 
event in the low momentum ir° region. In addition, the charged pion tends to be fast for the slow ir° region, thus 
leading to increased K* + <-> p + misidentification. 

In contrast, the only drawback of the fast ir° region over the three track sample is a factor two degraded AE 
resolution. We therefore choose to use only the half of the sample that has a high momentum tt° in our fits in the 
two track ir° final state at this point. Besides this, the same likelihood fits are made as described for the three track 
final state. 

Efficiencies and results are summarized in Table IV. The crossfeed rates from other PV modes as well as b — > c 
decay backgrounds are negligible. A significant signal in B° — > 7r ± /9 =F is observed at a branching fraction of B(B° — > 
tt^p^) — (3.5 ± 0.5) x 10~ 5 . Note that we do not tag the flavor of the B in the present analysis. The measured 
branching fraction is therefore the sum of B° — ► p + n~ and B° — ► p~n + . In addition, averaging over charge conjugate 
states is as always implied. 



Contour Plot 




10 20 30 40 50 

B°^-7T + p" Yield 



p — > 7T 7T° Mass 



FIG. 5. Contour and projection plot B° — » n ± p^ 1 



3. Evidence for B° -> tt~ K* 



We search for B° — > ir~K* + with submodes K* + — ► K^tt + and K* + — > K + ir°. Due to the large combinatoric 
and physics backgrounds in the soft ir° region, we only select the hard ir° region for the K + ir° decay of the K*. 
Backgrounds other than those from continuum are negligible. Event selections are presented in Table II. Efficiencies 
and results are summarized in Table IV. 

The individual branching ratios obtained in the two K* + submodes are consistent, and we combine the two submodes 
to arrive at an average branching ratio of B(B° — > n~K* + ) = (2.2 g in's) x 10~ 5 which is 5.9<r from zero. We note 
that the statistical significance depends largely on the two track K® final state, which has less background and larger 
efficiency than the two track ir° final state. In contrast to the two observed B — * pir decays, the one dimensional 
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projections of the fit (see Fig. |6[) are somewhat less than inspiring, and a simple event count in the mass plot does 
result in an excess of only 2.4<r. However, goodness of fit (21%CL), and likelihood per event distributions are perfectly 
consistent with expectations from Monte Carlo. The most likely signal events have signal likelihoods consistent with 
what one may expect from signal Monte Carlo, rather than the background data taken below BB threshold. 

In addition, we generated 25000 distinct Monte Carlo background samples in the K^tt + tt^ final state. Each of 
these samples has the same number of events as our actual data in this final state. We perform a likelihood fit to each 
of these 25000 samples and record signal yield and significance as reported by each fit. We find that none of these 
background samples leads to a reported yield or significance as large as found in data. We therefore conclude that 
our result is exceedingly unlikely to be due to a background fluctuation. 




TABLE IV. Summary of CLEO results for B decays to a pseudo-scalar and a vector mesons (PV modes) 





Mode 


Eff (%) 


Yield 


Signif 


BR/UL (10" 5 ) 


n ± p il 


30 ±3 


26.ll y 8 . 


5.2a 


1.5 ±0.5 ±0.4 




12 ± 1 


28.5 +8 7 9 9 


5.6a 


3.5l {'I ± 0.5 




7±1 


10.81" 


5.2a 


2.3±g : | ±0.3 




4.1 ± 0.4 


5.7±*;i 


2.5a 


O A+ 1 '' 5 + - 3 
z - u -l.l -0.4 






5.9a 


9 9+O.8 +0.4 
z - z -0.6 -0.5 




18 ±2 


12.31$;? 


< 3a 


< 2.7 


K±p° 


28 ±2 




< 3a 


< 2.2 


K°p° 


10 ± 1 


o 0+4.9 
°- Z -3.9 


< 3a 


< 2.7 


K±pT 


11 ± 1 


8.3±|;g 


< 3a 


< 2.5 




26 ±3 






< 0.59 


K°cf> 


7± 1 






< 2.8 




26 ±3 






< 0.40 




17 ± 2 






< 0.54 




7± 1 


o o +0 - 9 

u - u -o.o 




< 0.8 




4.1 ± 0.4 


o o +1 - 3 




< 1.7 


K ± K *T 








< 0.6 


K+K*°(K+7T-) 


18 ±2 


0+ 21 




< 1.2 
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IV. DISCUSSION OF OUR RESULTS 

Let us start by summarizing some of the more striking features seen in the data. First of all, we see no evidence 
for B — > KK decays in either B — ► PP or B — > PV . Such decays would proceed either via highly suppressed 
W— exchange (e.g. B — > K + K~) and & — > d penguin diagrams (e.g. _B — > K^K^, B — > K®K®) or via final state 
rescattering (FSI) . Given that our upper limits for some of these decays are an order of magnitude smaller than at 
least some of the branching fractions we measure it seems fair to neglect FSI when trying to understand the dominant 
contributions to charmless hadronic B decays. 

Second, we see no evidence for B — > tttt decays while we observe both B — > Kir as well as B — > pir decays. We try 



to make sense out of this in Section [V A in the context of isospin and factorization. 

Third, we are so far unable to measure the branching fraction for any of the B — > pK decay modes, despite the 
fact that we have measured B — > pir and B — > Ktt, and at least one of the B — > K*tt decay modes. This is in 
full agreement with factorization predictions. Factorization predicts destructive (constructive) interference between 
penguin operators of opposite chirality for B — > pK (B —* Kir), leading to a rather small (large) penguin contribution 
in these decays. In addition, factorization and CVC predict that only the left-handed penguin operator contributes 
in B — > K* + n~ . Destructive interference of penguin operators is therefore not expected in this decay mode. 

Fourth, we want to note that the measured ratio R p = B(B° — ► p ± n T )/B(B + — > p°tt + ) is much smaller than naively 
expected. In B — > pit decays the p can either come from the upper or lower vertex, and it is generally believed that 
upper vertex p production clearly dominates due to favorable form factors as well as decay constants. In addition, 
B + — ► p°ir + is further suppressed by a factor two because only the uu part of the p° wave function contributes. The 
present CLEO measurement of B(B° — > j o ± 7r T ) is the sum of upper and lower vertex p production. It is therefore 
rather surprising that the measured R p = 2.3 ± 1.3 is not significantly larger than two. Measurements of B — > p + Tr° as 
well as a flavor tagged measurement of B — ► p + ir~ would help to clarify the situation in B — ► pir decays. It remains 
to be seen whether or not such measurements are within reach using the full CLEO data set. 

Finally, maybe the most striking observation in our data are the large branching fractions measured for charged as 
well as neutral B decays to rfK . Violation of a sum-rule proposed by Lipkin jlj| seems to indicate that a significant 
flavor singlet contribution is needed to explain these rates. The literature is full p9] of attempts to explain this 
apparent discrepancy, the most interesting of which is the suggestion that R-parity violating couplings may explain 
the large rj'K as well as the stringent limit on (f>K J21J. The latter is particularly amusing as one of the relevant 
couplings (A323) would also be present in B s ~ mixing j2^| and could therefore lead to a different value for 7 as inferred 
from B — ► Kir decays and the limit on Am s / Amj in the context of the usual analysis of the p — r\ plane p0||. 



A. Understanding the non-observation of B — > nn 

Most theoretical predictions lead us to expect a branching fraction for B — * ir + ir~ at a level of 1 — 2 x 10~ 5 . pjj 
Instead, the central value and 90% confidence level upper limit presented here are 4 and 8 x 10 -6 . With results like 
this a natural question to ask is "How small can B(B — > 7r + 7r~) be?\ 

Let us start our answer by describing a data based factorization prediction. Assuming factorization, and neglecting 
W-exchange, penguin annihilation, and electroweak penguin diagrams one may expect the following expressions for 
the decay amplitudes: 

V2A ±Q = -(T + C) 
A+- = -(T + P) = -|T|e ,; T x (1 - \P/T\e ia ) (2) 
V2A 00 = P~C 

Superscripts +,—,0 indicate the charge of the final state pions, and T,C,P stand for external and internal W- 
emission, and gluonic penguin diagrams respectively (Fig. 1(a), Fig. 1(c), and Fig. 1(b)). 
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We can arrive at "data based factorization estimates" of these amplitudes if we identify C — 0,2/0,1 x T and use 



02/01 = 0.21 ±0.14 from measurements in B — > D decays 25 . We then estimate T using factorization and the CLEO 
measurement B(B -> nlv) = (1.8 ± 0.5) x 10" 4 @ as follows: 



T ~ Vfafir x at x W rCB^Tri^) x yJB{B^-Klv) 

~ l.OGcV x (1.0 ±0.1) x (0.27±0.05)/GeV x (0.0135 ± 0.0022) ^ 
- (3.6 ±0.9) x 10~ 3 

The dominant error here is due to the spread among a variety of theoretical models for the q 2 dependence of the 
form factor j27j. We do not assign any error due to a possible breakdown of the factorization hypothesis. Throughout 
this paper we express the absolute size of amplitudes in units of ^Branching Fraction. 

The decay B + — ► K® tt + has three down type quarks in the final state. Inspection of Figure 1 shows that this final 
state can only be reached via penguin diagrams, or final state rescattering. Furthermore, the electroweak penguin 
contribution to this decay is color suppressed, rather than the color allowed one shown in Figure 1(d). It is therefore 
reasonable to estimate P from the measured B(B — > K ^) corrected by CKM and SU(3) breaking factors. 

Using these numbers we arrive at \T/P\d = 5.0 ± 2.3. This leads to the factorization predictions B(B° — ► ir + ir~) = 
(8 ± 5) x 10~ 6 , B(B+ -> tt+tt ) = (10 ± 5) x 10~ 6 , and B(B° -> 7r°7r°) - few xl0~ 6 . The last of these three 
estimates is not very meaningful given the errors on the quantities that enter. We assume maximum destructive 
interference (cosa = 1). Ignoring the penguin contribution (i.e. cosa = 0) leads to a prediction of B(B° — > ir + n~) = 
(13.0 ±6.5) x 10~ 6 . 

As an aside, we can calculate \T/P\ S = 0.26 ± 0.08. This means that CP violating rate asymmetries as large as 
50% are in principle possible for decays like B — > K + ir~ if the relevant weak and strong phases are close to ±tt/2. 
In addition to these factorization estimates, it is quite illustrative to look at the isospin decomposition of B — > irir: 0] 

V2A ±0 = 3Aj /2 x e lS 

A \ Q = A 3/* Xei5 s +A l/2 + A fy W 
V^A 00 = 2A^ 2 x e tS - A\i 2 - Ay 2 

Here the subscripts 1/2,3/2 indicate the two different isospin amplitudes. Note that only the A 1 / 2 amplitude has 
any contribution from b — > d penguins, whereas the A 3 / 2 amplitude is a pure b — > u transition. We indicate this by 
making the dependence on weak (/3, 7) and strong interaction phases (5) explicit. Using the factorization estimates 
above, it is easy to show that |A^ 2 |, |-Ag/ 2 |, and \A^ 2 \ are of the same order of magnitude. 

Equation ^ shows that B(B — > tt^tt ) can be estimated without making any assumptions about strong or weak 
phases. However very little can be said about the relative size of B — > 7r°7r° versus B — > 7r + 7r~ without making such 
assumptions about relative phases. Common prejudice assumes 5 << 1 and therefore B — > ir + n~ >> B — > 7r°7r° due 
to the destructive interference between and Ai/ 2 in B — > 7r°7r . However, as we allow for S to increase towards it 
we not only decrease (increase) B — > tt + it~ (7r°7r°) but also increase the size of the "penguin pollution" in any future 
attempt of measuring sin 2a via time dependent CP violation in B — > ir + n~ . 

We are thus in the amusing situation that we would like B — > 7r°7r° to be large to make the Gronau, London isospin 
decomposition experimentally feasible. Though at the same time, we can only hope for 5 << 1 (i.e. vanishingly 
small B — > 7r°7r°) to avoid destructive interference between the two b — > u pieces in the amplitude for B — > ir + Tr~ . 

We conclude that our present data is still consistent with factorization predictions for B — ► 7r + 7r~. However, 
B — ► 7r + 7r~ could be significantly smaller than predicted by factorization if the strong interaction phase between 
isospin amplitudes is non-zero. 



*We ignore a possible strong phase difference between penguin and tree contribution to A 1( / 2 . 



12 



B. Comment on Neubert-Rosner bound on 7 



As previously mentioned, the ratio = B{B ± — > K°ir ± )/2B(B ± — ► K ± tt°) ||, may be used to constrain COS7 
if R± 1. Our measurement of this ratio is i?* = 0.47 ± 0.24. The relevant equation for bounding COS7 out of the 
paper by Neubert and Rosner || is: 

C0S7 = Sew - ((1 - \/#*)/ e 3/2)/ cos cj> + 0{e\ /2 ) (5) 

The parameter e 3 / 2 is defined in terms of experimentally measurable quantities below. It is essentially given by the 
ratio of b — > u tree and b —> s gluonic penguin amplitudes. The 0(e^ 2 ) terms were shown to be small in Ref. [ p8| . 
Here, Sew — 0.63 ± 0.15 is the theoretically calculated contribution from electroweak penguin operators ||. 

The dominant uncertainty in Equation ^ is the unknown strong phase cos 4>. Taking the extreme values of and 7r 
for this phase we thus arrive at an excluded region for cos 7 rather than an actual measurement. To be conservative, 
one may choose values for Sew such as to minimize this excluded region: 

0.48 + |1 - > cos 7 > 0.78 - |1 - y^K\/e 3/2 (6) 

The structure of this is obviously to exclude values for COS7 near 00S7 = Sew if -X" = |1 — v / ^*l/ e 3/2 > 0.15. The 
size of the exclusion region is determined by the central value of X as well as its error. The variable X defined here 
is given in terms of measurable quantities up to small uncertainties due to non-factorizable SU(3) breaking: 

|1 — a/b\ x a/c 
y/B(B -y K°n+) 
(3.74 ±0.72) x 10~ 3 

^2 x BjB -> K+tt°) (7) 
(5.48 ±0.91) x 10~ 3 
\V us /V ud \f K /U x ^2 x BjB -+ n+ir ) 
(0.95 ± 0.31 (0.91 ± 0.18)) x lO" 3 

The two different values for c are obtained using either the most likely value for B(B — > 7r + 7r°) based on the 
preliminary CLEO results (including statistical and systematic errors added in quadrature), or a weighted average 
of the latter with theoretical predictions based on factorization ||. When calculating X from these numbers we 
additionaly increase c to conservatively account for theoretical uncertainties due to non-factorizable SU(3) breaking 
("/if//4' = 1-33 rather than the experimental value of 1.2). The resulting values for X are 1.15±0.62 and 1.20±0.56 
respectively for the two different values for c. In the following we will use X = 1.20 ± 0.56. A number of comments 
are in order at this point. 

First, the central value for X leads to a physical value for cos 7 via Equation || only if the strong phase 0^0. In 
that case, the measured X then prefers rather large values of I7I ~ 120°. Such values of 7 are generally not the favored 
ones as they would tend to imply B s mixing to be smaller than the present limits and/or fg a y 7 Bb 3 / fB d y 1 ^B d to 
be at the large end of the generally assumed range. It was pointed out by He, Hou, and Yang |2j| that a number of 
other charmless hadronic B decay results from CLEO also suggest I7I > 90°. 

Second, only 10 — 15% of a Gaussian with mean 0.48 + X and u = <Jx ly within the physically allowed region for 
cos 7. Calculating a bound in this case isn't all that meaningful. Instead one may consider cos cf) < to be ruled out 
at ~ 90% confidence level. Using the usual procedure of calculating one-sided confidence levels based on the area 
inside the physical region only, results in the bound cos 7 < 0.33 @ 90% confidence level. 

Third, the experimental errors on X are large, roughly 1/4 of the physically allowed region total. It is fair to say 
that the only reason why we may deduce a non-zero exclusion region for cos 7 from present measurements is because 
our present central value for X indicates a prefered value for cos 7 that is far away from cos 7 = Sew- This is in 
contrast to some of the recent analyses of the p — rj plane J30| which tend to favor cos 7 ~ Sew- 



X=\l-^m\/e V2 ) = 
a =: 

b =: 
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V. CONCLUSION 



In summary, we have measured branching fractions for three of the four exclusive B — > Kit decays, as well as the 
two B — > rf K decays, while only upper limits could be established for all other B decays to two pseudo-scalar mesons. 
In addition, we have observed two of the four B — > pir decays, as well as one of the four B — > K*tt decays. We do not 
observe significant yields for B decays to pK, K*K , cfm or <j>K. 

The pattern of observed decays is broadly consistent with expectations from factorization. We see significant 
contributions from both b — ► u as well as b — + s transitions. 

In addition, the Neubert-Rosner bound derived from present CLEO data on charmless hadronic B decays indicates 
cos 7 < 0.33 @ 90% confidence level. This is in slight disagreement with some of the more aggressive analyses of the 
p — r\ plane found in the literature which prefer larger values of cos 7. 

Many thanks to our colleagues at CLEO for many stimulating discussions as well as the experimental work that 
made this paper possible. Further thanks go to A. Ali, J.-M. Gerard, M. Gronau, W.-S. Hou, M. Neubert, J. L. Rosner, 
and H. Yamamoto for discussions on topics related to Section IV. We gratefully acknowledge the effort of the CESR 
staff in providing us with excellent luminosity and running conditions. 
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